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Abstract 



We demonstrate the feasibility and efficiency of the Coulomb-Sturmian 
separable expansion method for generating accurate solutions of the Fad- 
deev equations. Results obtained with this method are reported for sev- 
eral benchmark cases of bosonic and fermionic three-body systems. Correct 
bound-state results in agreement with the ones established in the literature 
are achieved for short-range interactions. We outline the formalism for the 
treatment of three-body Coulomb systems and present a bound-state calcula- 
tion for a three-boson system interacting via Coulomb plus short-range forces. 
The corresponding result is in good agreement with the answer from a recent 
stochastic- variational-method calculation. 
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I. INTRODUCTION 



Separable expansion schemes have always been extremely useful in solving few-body 
problems. For example, in the three-body systems, at any time starting from the 60's, 
important new results have been achieved by solving the Faddeev equations with some sort 
of separable representation of the two-body subsystems (for a review see, e.g., ref. In 
more recent times, above all the separable expansion method proposed by Ernst, Shakin, 
and Thaler |0] (the so-called EST method) has proven very useful. Indeed, the first three- 
nucleon scattering results with such realistic meson-exchange N-N interactions like the Paris 
potential were achieved along this method |^ , and they were later on confirmed by a direct 
solution of the Faddeev equations @]. 

Even nowadays when three-body Faddeev equations can directly be solved on supercom- 
puters (for a recent review of the state of the art see ref. 0]), separable expansion schemes 
have their relevance. Not only does an accurate separable representation of the input dynam- 
ics allow to save much computer time in arriving at standard results, separable expansion 
methods may also help a lot in obtaining solutions to hitherto unsolved problems. In this 
respect we may mention the solution of the three-nucleon scattering problem with realistic 
N-N interactions and Coulomb forces at any energy; so far only limited solutions below or 
above breakup threshold have been obtained (see, e.g., refs. 0] and 0]). 

In this paper we deal with a separable expansion method that is well-adapted to treating 
few-body problems including long-range forces. Its essence lies in the expansion of the 
potential operator of the short-range part of any interaction with the use of Coulomb- 
Sturmian (CS) functions If the full potential also contains a Coulomb interaction , 
this is kept in the Green's operator. Thereby all difficulties associated with a (separable) 
expansion of the Coulomb potential are avoided, while at the same time correct asymptotic 
properties of all quantities are guaranteed. Still the advantages of the separable expansion 
can be exploited in solving the two- and three-body integral equations. 

The CS separable expansion method has been extensively tested before in two-body 
problems. Not only bound and resonant states with a variety of short-range plus Coulomb 
potentials were investigated but also scattering solutions were obtained . In these works 
also convergence studies were performed and subsequently extended to the multichannel 
Coulomb problem [|TT|]. Computer codes for the CS separable expansion of any local or non- 
local two-body interaction under the presence of Coulomb-like potentials were published 
in ref. ||lO|. In Ch. II below we shall recall some of the most important formulae for the 
two-body problem. 

The principal advantage of the CS separable expansion is the fact that the matrix el- 
ements of the Coulomb Green's operator can be calculated analytically in the two-body 
system. Besides the separable representation of the short-range part of the full interaction 
this turns out to be an essential requirement for an efficient and accurate solution of the 
three-body system. For the latter we may thus follow the integral-equation approach and 
thereby guarantee for the implementation of the appropriate (Coulomb- like) asymptotics. 
In Ch. Ill below we shall demonstrate how the matrix elements of the three-body Coulomb 
Green's operator can be calculated in a reliable way. 

We prove the efficiency of our method through the solution of the Faddeev equations for 
three-body bound states interacting via various short-range forces and a case with additional 
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Coulomb interaction among all three particles. We adhere to problems for which benchmark 
results from other methods have been obtained already. It is found that in all cases excellent 
agreement is achieved. The method therefore appears promising as an efficient tool for 
solving three-body systems, as it can be adapted to more general cases (including three- 
body forces) and extended to scattering under the presence of long-range interactions. 



II. COULOMB-STURMIAN SEPARABLE EXPANSION 

We give a short account of the formalism of the CS separable expansion in the two- 
body system with short-range plus Coulomb interactions. Here we specify our notation and 
provide the formulae that are later on needed in the solution of the three-body problem. 



A. Basis functions 

The short-range potential operator in some angular-momentum state / will be ex- 
panded on the basis of Coulomb-Sturmian functions 



{r\nl) 



(n + 2/ + l)! 



1/2 



{2bry+^e-'"^L^J+\2br), 
(n = 0,l,2,...), 



(1) 



which are the solutions of the Sturm-Liouville problem of the hydrogenic system [jT2[. Here, 
2^2«+i rgpj^gggnt the Laguerre polynomials and b relates to the energy in the Sturm-Liouville 
equation. We take 6 as a fixed parameter thus working with energy-independent CS func- 
tions. They form a complete set 



N 

1 = lim y2 \nl){nl\ = lim l^v, 

n=0 



(2) 



where 



{r\nl) = -{r\nl). 



(3) 



With the N-th order unit operator In in Eq. we can now expand the short-range 
potential operator in the form 



N 



lim InViIn = lim \nl){nl\vf\n'l){n'l\. 



(4) 



n,n'=0 



If remains finite, we end up with a rank-N separable approximation. As a consequence 
the two-body problem can then be solved by algebraic methods [Q. 
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B. Short-range plus Coulomb interactions 



Let us now assume a two-potential case of short-range plus Coulomb-like interactions 

V = v'' + v'^ (5) 

and consider the homogeneous Lippmann-Schwinger equation for the bound state | ■?/'/) in 
some partial wave / 

\i^i)=9?\Emi). (6) 

Here gf{E) is the two-body Coulomb Green's operator 

gfiE) = {E - - v^r' (7) 

with the free Hamiltonian denoted by /i^. Using the expansion in Eq. (|]) one arrives at 
a linear system of homogeneous quations for the wave-function coefficients Ai^ = {nlltpi): 

[i9fm-'-mi = o. (8) 

It has a unique solution if and only if 

det[{gf{E))-'-vn=0. (9) 
The matrices involved are made up from the elements 

gZjE) = {nl\gr{E)\^l) (10) 

and 

vl^, = {nllvtlri'l). (11) 

While the latter matrix elements may be evaluated (numerically) for any given short-range 
potential either in configuration or in momentum space, the matrix elements of the Coulomb 
Green's operator between CS states can be calculated analytically 0; the corresponding 
computer code is available from ref. |T^. This fact then also allows to calculate the matrix 
elements of the full Green's operator in the whole complex plane, 

^{E) = {{gf{E)r'-v]r\ (12) 

what will be needed later on in the solution of the three-body problem with charged particles. 
After solving Eq. (H) for the coefficients the bound state lipi) can be expressed as 

N _ 

\i^i) = Y.Bjn9?'mnl), (13) 

n=0 

where the new coefficients result from the matrix multiplication B^i = v_iAi- We note that 
expression (|13D is distinct from a usual expansion of the state \ipi) with certain test functions. 
The explicit occurence of the Coulomb Green's operator always ensures the correct asymp- 
totic behaviour 0. This is an immediate consequence of the fact that only the short-range 
potential (but not the wave function) is expanded. 
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III. SOLUTION OF THE THREE-BODY BOUND-STATE PROBLEM 



We now extend the CS basis to the three-body system and demonstrate the solution of 
the Faddeev equations for bound states of three particles with any short-range interactions 
and under the presence of Coulomb forces. 

A. Short-range interactions 

The integral equations for the three Faddeev components \E'q, of the bound-state wave 
function read: 

|vi/„) = G„(i?)K|vi/^)+<|vi/,)] (14) 
with a,/?,7 a cyclic permutation. Here the channel Green's operators are defined by 

G^iE) = iE-H'-vX\ (15) 

where is the free three-particle Hamiltonian and the short-range interaction of the 
pair (/3,7). In the angular-momentum representation (omitting the explicit spin and isospin 
dependence from our notation) we define the CS basis for the expansion of the short-range 
interactions in the three-particle system as 

\ni'l\)a = \nl)a ^ \l^X)a, ('^5^^ = 0,1,2,...), (16) 

with the CS states from Eq. (|l|). Here / and A denote the angular momenta of the two-body 
pair (/3,7) and of the third particle a relative to the centre of mass of this pair, respectively. 
In the three-particle Hilbert space we have 

N 

1 = J™ a{nvl\\ = lim lN,a (17) 

n,u=0 

where the configuration-space representation in terms of Jacobi coordinates and rja reads: 

{^aVa\nulX)a = -^{CaVa\niylX)a- (18) 

After the CS expansion of the potentials v^, w^, and in the three-particle space, the 
Faddeev equations can be rewritten as 

1^7)]- (19) 

By applying the CS states Q,(ni//A| from the left, Eqs. (|I9|) turn into a linear system of ho- 
mogeneous equations for the coefficients of the Faddeev components Ai^^^^^^ = o(^^z^/A|\E'a): 

[{G{E))-'-v]A = 0- (20) 
A unique solution threreof exists if and only if 
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det[{G{E))-^ - v] = 0. (21) 
The matrices G{E) and v have a block structure and the matrix elements are given by 

m^\^nu,l'^X'^n'u' = (1 - Safs) a {uiyl X\v'Jn' u' I' X' ) p (22) 

and 

G^^x^n.,i'^x'^n'u'{E) = ^{niA\\G^{E)\n'^'\')^, (23) 

respectively. Notice that the matrix elements of the Green's operator are needed only 
between the same partition a whereas the matrix elements of the potentials occur only 
between different partitions a and (3. The latter may again be evaluated numerically either 
in configuration or momentum space. We have adopted the configuration-space version of 
the Balian-Brezin method pr5 |. 

For the calculation of the matrix elements of the Green's operator in Eq. (p3|) we proceed 
in the following way. We split the three-particle free Hamiltonian into 

H' = hl+hl, (24) 

i.e. the free motions in the two Jacobi coordinates. Then we define the two-body Hamilto- 
nian he as 

hu = hl+vl (25) 
Since the commutator of the Hamiltonians he and h^„ vanishes, 

set //q, / 

'h^,hl]=Q, (26) 
we may apply the convolution theorem by Bianchi and Favella in 



G^{E) = {E- h,^ - hlr' = ^J^de{E - e - h,J'\e - h^J-\ (27) 

Here the contour G encircles the spectrum of h^^ without penetrating into the spectrum 
of h^^ (cf. Fig. |l|). We note that in this integral the roles of h^^ and h^^ may also be 
interchanged. 

After sandwiching the above Green's operator between the CS states, the integral in Eq. 
(^) appears in the form 

Gi^x^nu,i'^x'^n'u'iE) = 7^ <f de ^{nl\{E - e - h^J~'\^')^ ^{iJx\{e - hl)-'\i7X')^, (28) 

^ 'TCt J C 

where the separate matrix elements occurring in the integrand are known from the two- 
particle case of the previous chapter (cf. Eq. (p!2D). 

After solving Eq. ( ^OD for the coefficients Ai^Xc^ny the Faddeev components can be ex- 
pressed as 

N 

|*a)= E B^n.G^{E)\nul\)^, (29) 

n,v=0 

where the new coefficients again result from the matrix multiplication B_ = v_A- As before 
we hint to the advantage that the representation ( |29| ) of guarantees for the correct 
asymptotic behaviour of the Faddeev component, due to the explicit occurrence of the 
Green's operator. This is of particular importance in the Coulomb case below. 
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B. Coulomb-like interactions 



In this section we extend the formulation of the three-body problem to the case of long- 
range interactions. We assume the subsystem interaction to be a sum of short-range plus 
(repulsive) Coulomb interactions 

= < + , (30) 



and adhere to the Faddeev equations in the form as modified by Noble |15| : 

|vI/,) = G^(E)[<|vI/^)+<|vI/,)] (31) 

with the Coulomb-like Green's operator 

G^JE) = {E~H'-v:-v^- - v^)-\ (32) 

Herein all long-range interactions for all subsystems are collected, in complete analogy to 
the two-body case (cf. Eq. (|)). 

As in the previous section the short-range potentials are expanded on the CS basis leading 
to an equation similar to (|l^). A unique solution thereof exists if and only if 

det[iG^{E))-' -v]=0, (33) 

where the matrices have the same block structure as before in Eq. (|2T|). The important 
point is that v contains only matrix elements of the short-range interactions; in fact they 
are completely equivalent to Eq. ([2^). 

For the calculation of the matrix elements of the Coulomb-like Green's operator we 
proceed along the lines of the two-potential formalism |T^. First we rewrite it using the 
resolvent equation 

G^AE) = G^E) + GAE){v'^ + < - n^)G^(i?) (34) 

with defined by 

G^E) = {E - H'^ - v: - - u''j-\ (35) 

Here we have introduced the auxiliary potential m^, which is required to have the asymptotic 
form 

C Za{Zp + Z^) 

Ua (36) 

as rja —>■ oo. It may be viewed as the effective Coulomb potential between the centre of mass 
of the subsystem a (with charge -|- Z^) and the third particle (with charge Zq,). The 
important role of the potential is that asymptotically it compensates the Coulomb tail of 
the long-range potentials + v^' in Eq. (Q). Thus the combination U = + v!!^ — u'^ can 
be subject to a separable expansion and effectively be treated as a short-range potential. 
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With the help of the formal solution of Eq. (0) we may now express the inverse matrix 

(G^^iE))-^ as 



(G^AE))-' = {G':{E))-' ~ U, (37) 

where U_ is constructed from the matrix elements 

Ui^x^nu,i'^x'^n'u' = a{niylX\{v^ + v^^ - M^)|nV'/'A% (38) 

resulting from the CS expansion. The matrix elements of the Green's operator G^[E) can 
again be calculated by the contour integral as before. We are then left with the integral 

^i^x^n.AKn'AE) = 7^J de M\{E - e - h^J-'\^')^ ^{i^X\ie - /iJJ-^|z7A')„, (39) 

^ Tit' >/ C 

where h^^ = /i^^ - and h^^ = /i^^ -v'^-v^. 

The calculation of the Faddeev amplitude |\1/q,) in Eq. (pTD is completely analogous to 
the short-range case of the previous section. Only in Eq. (|29| ) the Green's operator Ga{E) 
must now be replaced by the Coulomb Green's operator G^{E) of Eq. (|3^) . 



IV. TESTS OF THE METHOD 

In this chapter we demonstrate the performance of the method in calculations of various 
three-body bound states. We have selected cases for which benchmark results are already 
available in the literature. The comparisons will prove the efficiency of our method, especially 
in the situation when Coulomb forces are present. 



A. Illustration of the convergence of the CS expansion 

Before presenting the final results, let us demonstrate the convergence of the results 
for the three-body bound-state energies. For this purpose we take the example of the Ali- 
Bodmer (AB) potential fl^ between a-particles of mass M 



v'{r) = 500exp(-(0.7r)2) - 130 exp(-(0.475r)2) (40) 

without and 

v{r) = v%r) + 4e^/r (41) 

with Coulomb interaction. We use units such that h'^/M = 10.36675 MeV fm^ and e^ = 
1.44 MeV fm. 

Evidently the quality of the results will depend on the number of terms employed in 
the separable expansion of the (short-range) potential. We quote the values of the binding 
energies from calculations without (Table |) and with (Table |I|) Coulomb forces, taking into 
account different numbers of channels rich = 1,2,3,4 (corresponding to angular momentum 
states up to / = A = 0, / = A = 2, / = A = 4, and / = A = 6 employed). In all 
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cases it is observed that convergence up to 5 significant digits is comfortably acliieved witli 
= 20 terms applied for n and v in the separable expansion. Remarkably the speed of the 
convergence is everywhere similar, irrespective of how many angular-momentum channels 
are included and whether or not Coulomb forces are present. We note especially for the 
Coulomb case that the satisfactory convergence stems from reliable separable expansions 
of the potentials f and f/, which - from the point of view of scattering theory - are both 
short-range potentials; the fall-off of U is much slower than of u*, though. 

In principle, the convergence may also depend on the (range) parameter h of the Coulomb- 
Sturmian functions. We found, however, that the dependence is weak in a relatively large 
interval of possible choices, just as it was established in the two-body case 0,^]. 

The convergence is practically of the same quality in the case of the other potentials 
considered below. 



B. Results for various three-body bound states 



We now present our converged results for the Malfliet-Tjon (MT) [|T^ and AB |17 
potentials and compare them to other benchmark calculations. For the MT potential 

v\r) = Kexp(-3.11r)/r - Kexp(-1.55r)/r (42) 

between two nucleons of mass m we consider two cases 
MTVa: K = 1458.0470 K = 578.0890 
MTVb: Vr = 1438.4812 K = 570.3316 

(here we use units such that fi^ /m = 41.47 MeV fm^). For MTVa we may compare to the 
results of the Los Alamos-Iowa and Groningen groups, both of which are obtained 



from a direct solution of the Faddeev equations in configuration space, and in addition to 
results from an ATMS P2| calculation, a Green's function Monte-Carlo (GFMC) calcula- 



tion p3|, an integro-differential-equation approach (IDEA) p4[, and a stochastic variational 



method (SVM) From Table |IT| it is evident that our method provides very accurate 

predictions for binding energies in all cases. This is true for the channel-by-channel com- 
parison with the Los Alamos- Iowa calculation and likewise for the comparison with the best 
results from the other works. With respect to the best results quoted in the lower rows of 
Table |T| we note that for the corresponding calculations the number of angular-momentum 
channels employed is either not definitely known or not specified in a scheme like ours. 

In order to demonstrate that in addition to binding energies our method also provides 
accurate three-body wave functions, we calculated the root-mean-square (rms) radius < 
^2 >i/2 Corresponding results are given in Table 11^ for the MTVb potential, in which case 



we can compare to the calculations of the Los Alamos-Iowa group . For both the binding 
energy and the rms radius the channel-by- channel comparison indicates perfect agreement. 

In case of the MTI-III potential for a system of three fermions, acting in singlet and triplet 
states, we may compare to the 2-channels calculation of the Los Alamos-Iowa group. For the 



MTI-III potential as parametrized in ref. ||T9[, we obtain a binding energy (converged result) 
of E = 8.5358 in comparison to = 8.536 calculated by the Los Alamos-Iowa collaboration. 

Finally we come to the comparison of the results for the binding energy of a system of 
three bosons interacting via the AB and AB plus Coulomb potential (Table |Vp. We may 
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compare our 4-channels result to calculations with the ATMS method (uncharged case 
only) and the SVM [^,^]. We do not know of any Faddeev results in this case. Again we 
show predictions for the binding energies and rms radii. All the values quoted in Table [V| 
show a convincing agreement of our results with the ones from the other approaches. We 
specially stress the agreement of the result for the case of rigorously including the Coulomb 
interaction with the rather reliable answer from the SVM. 



V. CONCLUSION 

We have suggested a separable expansion scheme, relying on Coulomb-Sturmian basis 
functions, for solving the three-body problem. The method is especially suited to the case 
when Coulomb-like interactions are present in one or all subsystems. It allows to solve the 
three-body integral equations by expanding only the short-range part of the interaction in a 
separable form while keeping the effect of the long-range part in an exact manner via a proper 
integral representation of the three-body Coulomb Green's operator. As a consequence the 
method has good convergence properties and can practically be made arbitrarily accurate 
by employing an increasing number of terms in the separable expansion. The usage of 
the Coulomb-Sturmian basis is essential to allow for the accurate evaluation of the matrix 
elements of the Coulomb Green's operator. 

Beyond the studies of the method in systems with two-body asymptotics conducted be- 
fore [P-pJ|, we have now demonstrated its convergence properties and efficiency in (bench- 
mark) calculations of the three-body bound-state problem without and with Coulomb inter- 
actions. In both cases the solution of the Faddeev equations show a rapid convergence, and, 
whenever a comparison is possible to existing results in the literature, correct predictions 
for the binding energies and wave functions are achieved. 

The method is capable of treating any kind of short-range interactions, even in the case 
when Coulomb-like forces are present. The solution of the three-body bound-state problem 
was carried out here. However, the method is also applicable for scattering problems. In 
this regard it has been proven useful already in the two-body case p|-[TI|. To solve the 



corresponding problem for a (charged) three-body system with the Faddeev equations some 
technical details in connection with the evaluation of the then occurring matrix elements 
still need to be worked out. 
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TABLES 



TABLE L Convergence of the binding energy of a three-boson system interacting via the 
AU-Bodmer potential, Eq. (40), with increasing basis for the separable expansion. N denotes the 



maximum number of basis states employed for n and u in Eq. (19). 



N 


1 


Number 
2 


of channels rich 

3 


4 


12 


4.13899 


5.11756 


5.17699 


5.17888 


13 


4.14092 


5.11911 


5.17862 


5.18057 


14 


4.14044 


5.11897 


5.17846 


5.18043 


15 


4.14046 


5 11898 


5 17850 


5.18047 


16 


4.14065 


5.11917 


5.17871 


5.18069 


17 


4.14069 


5.11926 


5.17880 


5.18079 


18 


4.14069 


5.11926 


5.17881 


5.18080 


19 


4.14069 


5.11926 


5.17881 


5.18080 


20 


4.14070 


5.11927 


5.17882 


5.18081 


21 


4.14071 


5.11928 


5.17884 


5.18083 


22 


4.14071 


5.11929 


5.17884 


5.18083 


23 


4.14071 


5.11929 


5.17884 


5.18083 


24 


4.14071 


5.11929 


5.17884 


5.18083 





TABLE II. Same as in Table | 


] for the Ali-Bodmer plus 


Coulomb potential, Eq. (pT|). 


N 




Number of channels rich 






1 


2 


3 


4 


12 


1.90151 


2.81629 


2.86703 


2.86839 


13 


1.90473 


2.81833 


2.86912 


2.87054 


14 


1.90368 


2.81796 


2.86871 


2.87014 


15 


1.90373 


2.81797 


2.86875 


2.87019 


16 


1.90404 


2.81819 


2.86899 


2.87044 


17 


1.90397 


2.81824 


2.86904 


2.87049 


18 


1.90401 


2.81824 


2.86905 


2.87050 


19 


1.90400 


2.81824 


2.86905 


2.87050 


20 


1.90401 


2.81825 


2.86906 


2.87051 


21 


1.90402 


2.81826 


2.86907 


2.87053 


23 


1.90402 


2.81827 


2.86908 


2.87053 


23 


1.90402 


2.81827 


2.86908 


2.87053 


24 


1.90402 


2.81827 


2.86908 


2.87053 
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TABLE III. Binding energies in the case of the MTVa potential for a system of three bosons. 







1 


Angular momentum channels 
2 3 


4 


This work 
Faddeev |0| 


8.04251 
8.0424 


8.22953 8.24978 
8.228 8.249 


8.25215 
8.251 


Faddeev 
ATMS 1 
GFMC 
IDEA 1 

SVM 1^ 


[21| 
22| 
23| 
24 1 
!51 




8.25273 
8.26(1) 
8.26(1) 

8.25 
8.2527 





TABLE IV. Binding energies and root-mean-square radii for the MTVb potential for a system 
of three bosons. 







1 


Angular momentum channels 
2 3 


4 


-Eb 


This work 
Faddeev @ 


7.5398 
7.540 


7.7147 7.7338 
7.714 7.733 


7.7361 
7.735 



< >i/2 This work 1.7265 1.7117 1.7098 1.7095 

Faddeev @ 1.727 1.711 1.710 1.710 



TABLE V. Binding energies and root-mean-square radii in the case of the Ali-Bodmer and 
Ali-Bodmer plus Coulomb potentials (Eqs. (]40| ) and (0), respectively) for a system of three bosons. 



Ali-Bodmer 








Without 


With v'^ 


—Eb 


This wc 
ATMS 
SVM 1 


)rk 
2| 


5.181 
5.18 
5.18 


2.871 
2.872 


< ^2 >l/2 


This wc 
ATMS 

SVM 1^ 


)r 

2: 

25 


k 

^1 

1 


2.434 
2.43 
2.43 


2.517 
2.517 
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FIGURES 

FIG. 1. Contour C for the integral in Eq. (^) in case of the three-body bound-state problem. 
To the right of the dotted vertical line lies the (continuous) spectrum of /i^^ , to its left the (discrete 
and continuous) spectrum of /i^^ . 
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